Prom tunneling to photoemission: correlating two spaces 



A. A. Kordyuk/'^ V. B. Zabolotnyy,^ D. S. Inosov,^ and S. V. Borisenko^ 

^Institute of Metal Physics of National Academy of Sciences of Ukraine, 03142 Kyiv, Ukraine 
^ Leibniz- Institut fur Festkorper- und Werkstoffforschung Dresden, 01069 Dresden, Germany 

(Dated: November 1, 2005) 

Correlating the data measured by tunneling and photoemission spectroscopies is a long-standing 
problem in condensed matter physics. The quasiparticle interference, recently discovered in high- Tc 
cuprates, reveals a possibility to solve this problem. Application of modern phase retrieval algorithms 
to Fourier transformed tunneling data allows to recover the distribution of the quasiparticle spectral 
weight in the reciprocal space of solids measured directly by photoemission. This opens a direct 
way to unify these two powerful techniques and may help to solve a number of problems related 
with space/time inhomogeneities predicted in strongly correlated electron systems. 
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The electronic inhomogeneity in high-Tc cuprates 
(HTSC) seen by scanning tunneling spectroscopy (STS) 
[l| has attracted much interest from the scientific com- 
munity because of its clear relation to a central prob- 
lem of HTSC-the mechanism of their evolution with hole 
doping from the antiferromagnetic insulator to supercon- 
ductor. Recent breakthrough in the development of STS 
technique to a level where the Fourier transformed (FT) 
STS images have been crystallized into well-defined sym- 
metric patterns has revealed the existence of regular 
inhomogeneities which could relate the HTSC problem to 
self-ordering phenomena in correlated electron systems. 
An important step in this direction have been made by 
identifying of the essential part of these inhomogeneities 
with rather trivial quasiparticle interference (QI) through 
an impurity scattering hypothesis [3j S which says 
that the FT STS function, FS'(r), is proportional to a 
joint density of states C(q) which, in turn, is just an 
autocorrelation (AC) of the quasiparticle spectral func- 
tion A(k). However, a reason for the remaining non-QI 
inhomogeneities, as well as the scattering hypothesis by 
itself, have yet to be understood. This is due, in part, to 
lack of deep understanding of the scattering phenomenon, 
but mainly due to an absence of a direct transition from 
the joint q-space of FT STS to the reciprocal k-space 
where the distribution of the quasiparticle spectral func- 
tion, A(k), is known from the angle resolved photoemis- 
sion spectroscopy (ARPES) with a tremendous accuracy 
0. The AC procedure gives an attractive possibility to 
compare the FT STS images to the AC ARPES data but 
evidently suffers from the absence of a direct transfor- 
mation from the joint to plain reciprocal spaces: Due to 
presence of experimental factors such as matrix elements 
the process of fitting ^(k) to C(q) is time consuming 
and, more important, does not give a unique solution 
0, 0] . Here we suggest a procedure to uniquely recover 
A(k) from S{v), shifting the problem of STS and ARPES 
correlation to ARPES domain. 

To understand the essence of the concept, let us rewrite 
the above mentioned impurity scattering hypothesis in a 



symbolic way: 

\FS\ ^ACA^C. (1) 
The AC is defined as 

AC A(k) = j ^(k)^(k + q)dk = C(q), (2) 

and is therefore an irreversible and time consuming op- 
eration. 

Nevertheless, we have found that both problems can 
be solved using the Wiener-Khinchin theorem Q, which 
sets an extremely important relationship between AC 
and FT: 

AC A = F-i|FA|2. (3) 

The fast Fourier routine essentially speeds up the calcula- 
tions, but more important is that the theorem shows the 
way to reverse the autocorrelation procedure and derive 
A{k) from S{y): With Eq. ([3]), the impurity scattering 
hypothesis can be rewritten as 

\FA\ ^ ^/WSI (4) 

and the problem is reduced to recovering of a function 
from its Fourier amplitude R{p) = \'R{p) \ = |FA(k)|. 

Fortunately, similar problem arose in applied optics 
and a number of phase retrieval algorithms (PRA) were 
developed |0] . These algorithms involve iterative Fourier 
transformation back and forth between the object and 
Fourier domains with application of the known con- 
straints. 

Fig. 1 (top row) shows A(k) — a gaped quasiparticle 
density of states of an optimally doped CuO-bilayer — as 
well as the real part of its Fourier image, Re(FA), and the 
result of its autocorrelation, AC A. A(k) is simulated at 
20 meV binding energy based on experimentally deter- 
mined parameters of the bare band dispersion for an op- 
timally doped Bi-2212 compound [13] taking into account 
a d-wave superconducting gap A = Ao(cosA:a; — cos ky)/ 2 
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FIG. 1: Representations of quasiparticles in different spaces of high-Tc cuprates: the spectral function A in the reciprocal space 
k, the real part of its Fourier image FA in a joint real space p, and a joint spectral weight C — AC A in a joint reciprocal space 
q. The reverse transformation C — > A is possible by means of PRA (see Movie SI). The bottom row illustrates the stability of 
PRA to experimental uncertainty — a noise has been added to |Fyl| (see also Movie S2). ^(k) is simulated at 20 meV binding 
energy based on experimentally determined parameters of the bare band dispersion for an optimally doped Bi-2212 compound 
[iOj l taking into account a d-wave superconducting gap A — Ao(cosfca; — cosfcj,)/2 with Ao = 40 meV; k- and q- images cover 
exactly the 1st Brillouin zone. 



with Ao = 40 meV. We note, that Im(FA) = due to 
even symmetry of the spectral function, A(k) — A(— k), 
so, the problem can be reduced to recovering of a real 
smooth function from its modulus. It seems that such a 
procedure can be realized, in principle, but should evi- 
dently suffer from finite resolution and unavoidable un- 
certainty of the experiment. 

In turn, the PRA are essentially discrete. Their 
uniqueness, in case of positive and spatially confined ob- 
ject, has been proved theoretically while the stability to 
noise and speed of convergence are the issues of continues 
development (oj. The recovering of ^(k) from |F^| by 
means of PRA is illustrated by Movie 1 [ll||. Here we 
used a modified "input-output" algorithm the nth 
iteration of which can be formulated as follows: 
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if Re(^;) > 0, 
if Re(^;j < 0, 



(8) 



where R{p) is a "source" function, and /? is a constant 
which we choose between 1 and 2 as a compromise be- 
tween speed of convergence and stability of the algorithm. 
As an initial guess we used a Gaussian distribution with 
random noise: ^o(k) = exp(|kp/i«^) -t- noise. 

Bottom row of Fig. 1 illustrates the robustness of PRA 
to experimental uncertainty: A'(k) is recovered from a 
noisy |F^| (the noise has been simulated by a random 
value from a Gaussian distribution such that the stan- 
dard deviation of an infinite number of such values would 
be 50% of the average value of \FA\). The work of RPA 
in this case is illustrated by Movie 2 [ll|. 

Finally, we discuss the existent attempts to compare 
the STS and ARPES data in the q-domain 0, |E li- 
lt has been shown that the intensity maps measured 
by ARPES, when autocorrelated, does not result in so 
distinct spots as observed in FT STS images [![. The 
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correspondence can be made better assuming better en- 
ergy resolution in ARPES Q , or k-dependent matrix ele- 
ments in STS Q- We believe that both effects should be 
taken into account together with the gaped and highly 
anisotropic quasiparticle self-energy, but, in this paper, 
we do not purpose to define a region of parameters which 
result in spot patterns seen in FT STS. Instead, we sug- 
gest the approach to recover A(k) from FT STS shifting 
this problem to the ARPES domain. Comparing A(k) 
directly measured by ARPES to one determined from 
STS with the suggested approach will, in particular, clar- 
ify such long standing problems as tunneling matrix ele- 
ments [l2| and inconsistency in the values of the quasi- 
particle self-energies determined from STS and ARPES 
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experiments on high- Tc cuprates 

In summary, we have proposed a procedure to derive 
the distribution of the quasipartical spectral weight in 
plain reciprocal space of a crystal from STS measure- 
ments and demonstrated its robustness to experimental 
uncertainty. When applied to STS data, this procedure 
will give an ultimate check for the validity of the impu- 
rity scattering hypothesis and, if valid, open a direct way 
to unify STS and ARPES techniques. 

The project is part of the Forschergruppe FOR538. 
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